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Inverse Laplace transform on the lattice spacing is introduced as a computational framework of the extrapo- 
lation of the strong coupling expansion to the scaling region. We apply the transform to the two-dimensional 
non-linear O(N) model at N > 3 and show that the approximation of the continuum limit of the susceptibility 
agrees with the existing theoretical and Monte Carlo data. 
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The lattice spacing a serves as the fundamental parame- 
ter in lattice formulation of physical models. In addition to 
the role of ultraviolet regulator, it controls almost all physi- 
cal quantities, since the bare coupling g depends on it due to 
the dimensional transmutation. Further it is a direct measure 
of how the model being far from or close to the continuum 
limit. As an attempt to explore the theoretical framework of 
approaching to the continuum limit from the view point of the 
expansion around a — oo, we have studied the delta expan- 
sion and Pade-Borel technique in recent works GU^D- 

We are now aware that the techniques used in (ll 01 can be 
focused to the inverse Laplace transform on the lattice spac- 
ing a, which may be extended to complex valued for the sake 
of the transformation. The purpose of the present article is 
to give the basis of the method and, as an illustration, revisit 
the two dimensional (2D) O(N) sigma model at N > 3 to 
demonstrate the successful approximation of the susceptibil- 
ity in the continuum limit. 

As a conventional manner, it is convenient and natural to 
use a in dimensionless combination with other dimensional 
parameter. The appropriate one would be the scale but 
it is hidden with a inside the bare coupling. In a sense, the 
mass gap or the correlation length is a concrete realization of 
the scale A^ . They can be computed in the strong coupling 
expansion and the result fixes the functional relation between 
them and g at large a. Thus, for example, the mass gap (or 
square in some models) combined with a to lose mass dimen- 
sion is an appropriate parameter replacing A^a. Let us denote 
such a variable as M. 

Now we introduce the inverse Laplace transform as follows: 
For a quantity Q(M) given as a function of M, we make com- 
plex extension of M to z G C by adding pure imaginary part 
to M. Then, consider the superposition over the imaginary 
part and define inverse Laplace transform for t € K by 
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The contour of integration is parallel to the imaginary axis. 
Laplace transform recovers the original function by 



Q(z) = z dtcxp(~tz)Q{t) := L[Q}. 



(2) 



As long as the lattice models under consideration are subject 
merely to second or higher order transitions, there is no singu- 
larity at positive real part of z axis. Here, we further assume 
that Q(z) has no singularity in the right-half plane of z. This 
ensures that the contour position specified by M is arbitrary 
as long as M > and M disappears in the result. The first 
singularity to be met with would be, when we move the con- 
tour to the left, found at z = if the singularity coming from 
in the measure is not removed. 
As in the Fourier transform, the scaling behavior of Q(M) 
is connected to the large t behavior of Q(t) (see (Q~|i or ©). In 
turn, small t behavior is related to the large a behavior of Q. 
Then, as is apparent from the following basic result, 
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the series of Q in t is usually an entire one and quite adequate 
to approximate the scaling behavior of Q(M). The physical 
interpretation of t is not clear to us. However, the study of 
Q(i) in the approximate computation of the continuum limit 
of Q(M) is much more tractable than Q itself. 

Now, we apply our scheme to 2D non-linear O(N) sigma 
model at N > 3 and illustrate how working in the t-space 
endowed with Pade approximation is effective in the study of 
continuum limit via strong coupling expansion. The standard 
action of the non-linear sigma model on square lattice reads 



n fl=l,2 
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where ei = (1, 0), e.2 — (0, 1) and /3 is the reciprocal of the 
naive bare coupling g, 



P = i/g- 



(5) 



The vector a = (<ti, 02, • • • , <Jn) is constrained to satisfy 
<7 2 = N at all sites. We define the basic parameter M by 
the zero momentum limit of two-point function of a (see |Q]]). 
It was obtainable from the susceptibility and second moment 
both of which were computed in fi up to 11 For our ap- 
proach, j3 must be rewritten in M. It is simply given by the 
inversion giving |Q]] 
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At small M near the continuum limit, we notice the pertur- 
bative result of correlation length £ J3-0], 



£ = Cf exp 
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62 = 



2tt7V 
1 



N(N - 2) 
1 



(27T/V) 2 ' 

(-0.0490- 0.01417V), 



J 7V 2 (7V-2) 2 

x (0.0444 + 0.02167V + 0.00457V 2 - 0.01297V 3 ). 

The mass square M is the momentum counter part of £ 2 and 
its dependence on /3 should be identical with but possi- 
bly accompanied by another multiplicative constant, say Cm- 
However, Monte Carlo data |§] showed that the difference is 
less than a percent at TV = 4. Since the two constants agree 
with each other in the large TV limit due to the Gaussian nature 
of two point function, the difference between C^ and Cm may 
actually be negligible for all TV > 3. So, we neglect the slight 
difference of Cm and C^ and simply substitute M -1 / 2 into £ 
in 0, which is valid as long as B is large enough. The con- 
stant C(. is specified only non-perturbatively. It was exactly 
computed by Hasenfratz, Maggiore and Niedermayer ^ as 



C e = 32~ 1 / 2 ' 1 
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By solving (0 for B, we have 



^log*+^log ( ilog.)+ (^), (9) 
47r7V 27r7V 2 log a; 



where 
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Since the complete expression up to the four loop level and its 
transform are lengthy, we refer them in yl . The expansion (0 
will be used when physical quantities are expressed in M at 
weak coupling region. 

Representation of bare coupling in t-space is adequate for 
the extrapolation of the strong coupling series (0 to the scal- 
ing region. The inverse Laplace transform is obtained by 



changing M n to t n /n\, giving 
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Using the transformed bare coupling (fTTT l. one can confirm at 
least around or above TV = 6, the continuum scaling of /3, 
the inverse Laplace transform of (0. Then, the mass gap in 
the continuum limit can be estimated by fitting the four loop 
fipert to ( fTTT l. The process and the result is exactly the same 
as that presented in Hence we suffice ourselves to say the 
conclusion that the scaling behavior was observed in j3(t) in 
series of t and the estimation of agreed with the result (0. 

Now, we present a new attempt to compute the susceptibil- 
ity x near the continuum limit via its strong coupling expan- 
sion. Though the problem was discussed in IToll according to 
the technique of Pade-Borel approximants, B was taken as the 
function of \- The proposal in this brief report is to take a uni- 
fied approach in the study of the continuum limit by choosing 
A^a or M as the fundamental parameter. Here, we regard \ 
as the function of M and consider inverse Laplace transform 
of x on complex M plane. 

The susceptibility \ is defined by 
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From perturbative renormalization group at TV > 3, \ near the 
continuum limit behaves as 



X ~ (3C X exp 

where 77 = 2bi /6q 
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0.0626TV), 



TV 2 (TV- 2) 2 

x (0.1316 + 0.0187TV - 0.0202TV 2 



0.0108TV 3 ). 



The constant C x in dT3l > is of non-perturbative and not analyt- 
ically known. We estimate C x in later. 

First of all, we express \ m terms of the mass square M in 
both regions at 8 3> 1 and B <C 1. At weak coupling, using 
the perturbative result (0 and (fT3l . we find 



_C x (TV-2) f\ogx\ -i/(N-2) r 2(ai - loglogx) a 2 + a 3 logloga; + a4(logloga;) 2 
X= 2^N X V^2~) V + (TV-2) 2 logx + (N - 2) 4 (log x) 2 

aj = -2(TV - 2)n(2Nb 2 - c 2 ) + log 2, 

a 2 = 8(TV- 2) 2 tt[A(6 2 + 6TVtt^ -4TVtt6 3 -47r6 2 C2) + 27rc 3 ] - 4(TV - 2)(-l + 2(TV - l)7r(2TV6 2 - c 2 ))log2 

+2(A-l)(log2) 2 , 
a 3 = 4(7V - 2)(-l + 2(TV - l)n(2Nb 2 - c 2 )) - 4(TV - 1) log 2, 
a 4 = 2(TV- 1). 



(14) 
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We need to compute inverse Laplace transform of 
ar _1 (log z~ l )~ a (a > and z stands the complex ex- 
tension of M), We note that there is a singularity at z = 1 
coming from (log z)~ a . This an artifact of perturbation but as 
a consequence we need another cut in the z plane at [1, oo). 
Thus we restrict the location of the contour in the contour 



integral 
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L- 1 [z-\logz- 1 )- a }= / ^-z-^logz- 1 )-, 
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to be included within the strip < < 1. Perturbative 
region corresponds to the region ( » 1 and the asymptotic 
expansion of the integral reads 



L-'iz- 1 (log z- 1 )-"} = t(logt + lE )- a (l + 



a a(a + 1) 
1 — 

logt + ^E (logt + 7 £ ) 2 



(l-C(2)/2) + ...). 
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The transform of terms of the type 2(log z) "(log log z) k (k = 1, 2) can be readily obtained from ( TT5b by shifting a — > a + e 
and then expanding the result in e. In this manner we arrive at 
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where 



T = tjCle"*. 



(17) 



At strong coupling, the susceptibility is obtainable in pow- 
ers of p. The result of Butera and Comi gives us that 
X = /3(1 + 4/3 + 12/3 2 + 72 +^ 2 2 Ar /? 3 + •■■)■ Then, substi- 
tuting ^ into the series, we obtain the series of \ expanded 
in l/M, 
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X ~ M ~ (N + 2)M 2 
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The inverse Laplace transform is readily obtained by changing 
M -n _^ t n/ nl The result reads 



X = t 



2t 2 



2\{N + 2) 



+ 



2(597V 2 + 364iV + 480)t 3 
3!(A^ + 2) 2 (A^ + 4) 



Xstr- 



(19) 



To study the scaling behavior of \ from the series ( fl~9l . 
the series expansion of \ is not sufficient for accurate results, 
though inverse Laplace transform substantially enlarged the 
effective region of the truncated series. It requires expansion 
up to quite large orders which no one performed yet. For- 
tunately the status of approximation task is improved by the 
use of Pade approximants. The Pade approximant method 
replaces the power series by the ratio of polynomials in t. 
The degrees m and n of the numerator and the denomina- 
tor must have sum which agrees with the truncation order of 
Xstr- Since, at large t, it is perturbatively known that x be- 
haves linearly in t with the logarithmic correction, we con- 
fine ourselves with the rational functions at diagonal and near- 
diagonal entries in the Pade table. 
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FIG. 1. Pole distribution of diagonal Xstr for N = 3 at orders [3/3], 
[4/4], ■ ■ ■ and [9/9]. The gray level of dots represents the order. 
Lightest ones correspond to poles at 6th and black ones at 18th. 



Denoting the rational approximant as [m/n], we have ex- 
amined the pole distributions of [n/n], [n + 1/n — 1], [n — 
1/n + 1] and [n + n/n] cases. We found that [n/n] and 
[n + 1/n — 1] approximants have no poles in the right-half 
plane of t^ 1 for all N > 3. Further, the pole distribution 
from 4th to 18th orders at N = 3 shows that the poles in the 
left-half plane move toward the origin as the order increases 
(see Fig.l). It is highly conceivable that at the infinite order, 
the set of infinite poles collapse to the origin from the left. 
This is a reflection that the transformed series in t is an entire 
function and the set of degenerate poles represents the pertur- 
bative singularity implied by logarithms in ( fTo*l ). On the other 
hand, [n — 1/n + 1] and [n + 1/n] approximants are plagued 
with poles in the right-half plane. Even those poles are merely 
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artifacts of Pade approximant, they badly influence the behav- 
ior of X[n-iln+i\ on me rea l t-axi&. We therefore discard the 
result from [n — 1/n + 1] approximants. 




FIG. 2. Xpert (dotted) and X[9/9] f° r N = 3. In plot of Xpert, we 
used for C x the estimated value C x = 0.012104. 




C x = 0.012104 achieves the best optimization. Only at this 
value, both of Xpert and X[g/g] and their derivatives agree at a 
point. In this way, we can estimate unknown non-perturbative 
constant C x for other values of N by using [m/n] = [9/9] 
and [10/8] forms of approximants. The result is summarized 
in Table I. 

For all values of N at which Butera and Comi made the- 
oretical estimation of C x , our results are in good agreement. 
There is also Monte Carlo data for N = 3,4 and 8. At N = 3, 
the two results, C x = 0.0146(10) (H and C x = 0.0130(5) 
HI were reported. At N = 4, C x = 0.0329(16) in [U3n 
andCV = 0.0383(10) in HI. At N = 8, C x = 0.1037(4) 
in HU] and C x = 0.1028(2) in fH]. Though small discrep- 
ancy remains, our estimates are consistent with these Monte 
Carlo data. The values of [9/9] and [10/8] entries are close 

TABLE I. Estimated values of the non-perturbative constant C x . On 
estimation we have used Pade-Borel approximants of j3 at order 18th 
for the cases [9/91 and [10/8]. Cbt means the estimation due to 
Butera and Comi Q]. 



N 


[9/9] 


[10/8] 


Cbt 


3 


0.012104 


0.012034 




4 


0.037018 


0.036971 


0.034 


5 


0.059836 


0.059809 


0.059 


6 


0.078058 


0.078042 


0.077 


7 


0.092404 


0.092395 




8 


0.103830 


0.103824 


0.1035 


9 


0.113083 


0.113079 




10 


0.120703 


0.120700 


0.1212 


20 


0.157114 


0.157113 





FIG. 3. Plots of X[9/9] ~~ Xpert at N = 3 for three values of C x , 
C x = 0.011, 0.012104 and 0.013. 

After these observations, we compare Pade approximants 
at [n/n] and [n + 1/n — 1] to the 4-loop scaling behavior 
Xpert by using the exact result of C<t @|. We have ascer- 
tained that Pade approximants approximately recover the scal- 
ing behavior for N > 3. For example at N = 3, from the 
plot of X[9/9] and Xpert in Fig. 2, we see the two functions 
are close with each other for C x = 0.012104 in the interval 
around t ~ 40. To clarify the effect of C x on the behav- 
iors, we have plotted in Fig. 3 the difference Xpert — X[9/9] at 
C x = 0.011, 0.012104, 0.013. It is clearly seen that the value 



with each other for all N and this fact indicates the stability 
of the approximation. As a whole, the level of agreement with 
existing results is improved than the result of [ 10]. 

As demonstrated by the examination presented so far, in- 
verse Laplace transform on the lattice spacing is found to be 
effective in the approximation of the continuum limit. Rather 
than the real space of a, the dual space is convenient to con- 
sider scaling properties. The inverse Laplace transform pre- 
serves the essential information at scaling region and creates 
an entire function in t. It is worth of exploring the approach 
as an alternative route to the continuum limit from the strong 
coupling expansion. 
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